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> ■ ABSTRACT 

l> 
00 

C3 ' We report the results of ROSAT and ASCA X-ray observations of the supernova 

remnant N157B (or 30 Dor B, SNR 0539-69.1) in the Large Magellanic Cloud. For 



. comparison, we also briefly describe the results on SNR 0540-69.3, the only confirmed 

' Crab-like remnant in the Cloud. The X-ray emission from N157B can be decomposed 



into a bright comet-shaped feature, superimposed on a diffuse emission region of 
a dimension ~ 20 pc. The flat and nearly featureless spectrum of the remnant is 



O 

i-^ ' distinctly different from those of young shell-like remnants, suggesting a predominantly 

^ ■ Crab-like nature of N157B. Characterized by a power law with an energy slope ~ 1.5, 

the spectrum of N157B above ~ 2 keV is, however, considerably steeper than that of 
SNR 0540-69.3, which has a slope of ~ 1.0. At lower energies, the spectrum of N157B 
presents marginal evidence for emission lines, which if real most likely arise in hot gas 
\ of the diffuse emission region. The hot gas has a characteristic thermal temperature of 

0.4-0.7 keV. No significant periodic signal is detected from N157B in the period range 
of 3 X 10-3 - 2000 s. The pulsed fraction is ^ 9% (99% confidence) in the 2 - 7 keV 
range. 

We discuss the nature of the individual X-ray components. In particular, we 
suggest that the synchrotron radiation of relativistic particles from a fast-moving 
(~ 10^ km s~^) pulsar explains the size, morphology, spectrum, and energetics of 
the comet-shaped X-ray feature. We infer the age of the remnant as ~ 5 x 10'^ yrs. 
The lack of radio polarization of the remnant may be due to Faraday dispersion by 
foreground H II gas. 



Subject headings: galaxies: individual (Large Magellanic Cloud) — supernova remnants 
— X-rays: ISM 
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1. 



Introduction 



Only three supernova remnants (SNRs) have been unambiguously confirmed to be Crab-like: 
SNR 0540-69.3, MSH 15-52, and the Crab nebula itself. These remnants show centrally peaked 
morphologies and flat emission spectra; the emission comes predominantly from synchrotron 
nebulae powered by embedded young pulsars. There are a number of SNRs that have similar 
morphologies and spectral characteristics, but contain no detected pulsars. N157B is such a 
remnant (Mills, Turtle, & Watkinson 1978; Clark et al. 1982; Mathcwson ct al. 1983; Wang k 
Helfand 1991; Chu et al. 1992 — CKSL92 hereafter; Dickel et al. 1994 — D94 hereafter). Like 
SNR 0540-69.3, N157B is in the Large MageUanic Cloud (LMC) and thus has a good distance 
determination. But, previous observations have shown no clear correspondence between the 
morphologies of N157B in radio, optical, and X-ray. The remnant's age of ^ 2.4 x 10^ yrs, inferred 
from kinematics of optical features (CKSL92) is also much greater (by a factor of about 10) than 
the confirmed Crab-like SNRs. No radio polarized emission has been detected from N157B, as 
might be expected for a synchrotron nebula (D94). In short, the exact nature of the remnant has 
been uncertain. 

Here we report new results on N157B, based primarily on observations from the ROSAT and 
ASCA X-ray Observatories. These observations provide much improved timing, imaging, and 
spectral resolution capabilities. For comparison, we also present the results of analyzing a 40 ks 
ASCA observation on SNR 0540-69.3. Wc describe the reduction of the X-ray observations in §2, 
and present our analysis and results in §3. Then in §4, we discuss the implications of our results, 
together with previous observations, and postulate scenarios to explain various components of 
N157B. We summarize our results and conclusions in §5. Appendix A describes an algorithm 
we used for extracting ASCA spectra with correction for the instrument point spread function 
(PSF). Throughout the paper, we adopt an LMC distance D ^ 47 kpc (Gould 1995), at which 1' 
corresponds to 14 pc. 



Table 1 summarizes parameters of the two ASCA observations used in this study. Both 
observations were acquired during ASCA^s Performance Verification phase, and the data were 
made available through the public archive. Since the ASCA has two types of instruments, which 
observe simultaneously, a pair of Solid-state Imaging Spectrometers (SIS) and a pair of Gas 
Imaging Spectrometers (CIS), the two ASCA observations yielded eight sets of data: four from 
the SISs and four from the GISs. The SIS has an energy coverage approximately in the 0.5-8 keV 
range, and the GIS in the 0.8-10 keV range (Tanaka, Inoue, & Holt 1994). 



2. 



Observations and Data Reduction 



2.1. 



ASCA Observations 



Data from the SIS detectors were coUected in 4-CCD FAINT and BRIGHT modes. The data 
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were processed by converting the FAINT mode data to BRIGHT mode equivalent data and by 
filtering the combined set through the standard cleaning criteria. CCD pixels with anomalous 
counts were identified and excluded using the CLEANSIS technique (Gotthelf 1993). Detailed 
procedures of these methods can be found in Day et al. (1995). Fig. 1 presents an exposure 
corrected, smoothed SIS image of the 30 Dor region created by combining screened data from 
both observations. 

The SIS data are the most useful in our spectral study of N157B, because of both the 
broad-band energy coverage and the good spectral resolution {5E/E ~ 0.02(5.9 keV/Ef -^). The 
spatial resolution of the data is, however, limited by the poor point spread function (PSF) of the 
ASCA telescope. Although the PSF has a sharp core with a FWHM of ~ 1', the broad wing of the 
PSF produces a 50% encircled radius of ~ lf5 (Jalota, Gotthelf, & Zoonematkermani 1992). The 
PSF also contains some energy dependence which is most notable above ~ 6 keV. The broad wing 
of the PSF prevents an independent measurement of the background in the region of the source 
uncontaminated by the source spectrum itself. 

We thus devised a method that simultaneously solves for both the source and background 
spectra (see Appendix A for details). We applied this method to determine the ASCA spectra 
for both N157B and SNR 0540-69.3. To minimize the uncertainties in the spectral background 
subtraction, we chose a small on-source aperture of radius l'.7, and a concentric background 
annulus with its inner and outer radii of 2^3 and 4^9. We further used only the southwest half 
of the background annulus around N157B to avoid the possible contamination by the 30 Dor 
nebula in the northeastern half of the annulus. Within the on-source aperture of this remnant, the 
background accounts for about 7% of the total observed intensity with a relatively soft spectrum. 
Judging from the background intensity variation around N157B, the uncertainty in the background 
should be less than ~ 50%. This uncertainty in the background subtraction is probably within the 
systematic error of the data and the analysis. 

For each SIS data set of N157B we generated an appropriate background spectral file to 
match the on-source file. These spectra were then summed for both CCD detectors and for both 
observations. Individual response files were summed to create an appropriate count weighted 
multi-CCD detector and multi-telescope response matrix. The summation was accomplished using 
the FTOOL MATHPHA to add ASCA spectral files, and the combined instrument /response files 
were merged using the FTOOL ADDRMFQ. The spectra were binned to contain at least 25 counts 
for model fits. Our final spectral analysis of N157B was based on this summed spectrum and 
response. We also conducted fits to the four sets of spectra before the summation, and found that 
the results were consistent with those from fits to the summed spectrum, which was certainly 
much easier to analyze and made multi-component model fits possible. 

The CIS data are the most suitable for timing analysis to search for a putative pulsar in 



^ Information of the FTOOLS software package is available in the ASCA data analysis guide (Day et al. 1995). 
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N157B. All the GIS data, acquired in the highest PH time resolution mode (10 bit), sets a limit 
on the timing accuracy during space-craft telemetry down-links of 0.3 milli-seconds for HIGH 
bit-rate. However, the calibration of the absolute time (ASCATIME), is only considered accurate 
to about 2.0 ms (Hirayama et al. 1996, ASCA news no. 4). A pulsar search which combined high 
and medium bit-rate mode data limits the overall accuracy to 5 ms; low-bit rate mode data was 
excluded from our pulsar search. We converted the arrival time of each count to the barycentric 
dynamical time, using the FTOOL task TIMECONV. We also extracted a GIS light-curve from 
the ASCA observation of PSR 0540-69.3 to be used as a control on our procedure. 

The SIS data are of limited use in timing analysis, The 4 CCD mode SIS observation of 
N157B used a 16 s integration (time between CCD read-out), setting a fundamental limit on 
searches for frequencies faster than 0.03125 Hz, the appropriate Nyquist frequency. 

2.2. ROSAT Observations 

We utilized two ROSAT observations from the ROSAT High Resolution Imager (RHRI) and 
one from the ROSAT Positional Sensitive Proportional Counter (PSPC). Both the PSPC and 
RHRI were sensitive to photons in the ~ 0.1 — 2 keV range. 

The RHRI observations with a PSF of ~ 6" (FWHM) are the best for studying the spatial 
properties of N157B. One of the RHRI observations {ROSAT Seq. No. rh600228) with an effective 
exposure of 30094 s was pointed at R136, ~ 7' away from N157B; the other (wh500036) of 4460 s 
was pointed directly at N157B. Each RHRI observation covers a field of ~ 34' diameter. A coadded 
intensity image of the two observations is presented in Fig. 2. The processing of the observations 
was described in an earlier article by Wang (1995), where the same two observations were used to 
detect two point-like X-ray sources in the 30 Dor core region and to identify them as Wolf-Rayet 
-I- black hole binaries. Briefly, small astrometric errors {^^ 3") were corrected through comparisons 
between optical and X-ray positions of identified X-ray sources in the region. Therefore, the 
absolute astrometric accuracy of the observations is better than ~ 3". 

The PSPC observation (rp500131), targeted at R136 (Chu et al. 1993), was used to offer a 
useful combination of moderate spatial resolution (FWHM ~ 0^5) and spectral resolving capability 
(about six independent bands in the 0.1-2 keV range). The spatial resolution of the PSPC, 
though not as good as that of the RHRI, was still significantly better than the ASCA instruments. 
We extracted a spectrum of N157B from the PSPC observation, using an aperture of l'.5 for 
the on-source count rate and an annulus between 2'-3'.5 radii for the background estimation. 
The aperture is large enough to contain essentially all the counts from the remnant. The total 
background-subtracted count rate is 0.24 counts s~^ in the PSPC 0.1-2 keV band. Similarly, we 
obtained the RHRI count rate of 0.083 counts s~^ in the same aperture. 

We found that the PSPC and RHRI observations provided no strong constraints on the 
timing properties of N157B. The observations were taken in various time segments spreading over 
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a long period of time — about a half year in the case of rh600228. 

3. Analysis and Results 

3.1. Timing Properties 

We conducted an extensive search for pulsed emission from N157B. To familiarize ourselves 
with various methods used, we applied them first to SNR 0540-69.3, which contains a known 
pulsar. It was easy to detect the pulsar's period of 0.05041843 s, even with one tenth of the SNR's 
data. Above 2 keV, the pulsed fraction is 26 ± 2%, almost independent of photon energy. But the 

fraction falls to 15 it 2% in the 1-2 kcV range and to 9 ± 5% below 1 keV, apparently because of 
soft diffuse emission from the remnant. No significant periodic signal, however, was detected from 
N157B, confirming the preliminary conclusion reached by Itoh et al. (1994). In the following, we 
concentrate on describing our timing analysis of N157B with the GIS data. 

We made an FFT analysis of the combined light curve of the two GIS detectors. This analysis 
is sensitive to the period range from 3 ms to about 2000 s. The lower boundary is chiefly due 
to the number of time bins allowed by our computer memory, and the upper limit to various 
observing intervals or gaps on timescales comparable to, or greater than, 2000 s. A single FFT 
of the light curve is the most sensitive, but is limited to the period range longer than 0.05 s. 
For shorter periods, we first divided the total observing time into up to 16 equal time intervals, 
and then applied FFT separately to the individual light curves. The derived power spectra were 
averaged to form a combined spectrum, where we found no significant peak. We even tried the 
incoherent harmonics addition (e.g., Lyne &; Graham-Smith 1990) to strengthen harmonics, in 
case the pulse shape was narrow. But no consistent harmonic pattern was discovered. 

We further utilized a period folding method, but again detected no significant signal. A 
fit with a sinusoidal waveform to a folded pulse profile at an FFT peak period typically yielded a 
pulsed fraction of 6%. By re-sampling counts according to the "best-fit" model, we estimated 
a 99% upper limit to the pulsed fraction as ~ 9% in the energy band 2 — 10 keV. The limit for a 
sharper pulse profile would be tighter. A similar analysis with the PSPC data gave an upper limit 
as ~ 17% in the 0.5-2 keV band. 

3.2. Spatial Properties 

Figs. 1 and 2 provide global views of the region including both N157B and the 30 Doradus 
nebula. The RHRI image shows that N157B is well separated from the main body of the 30 
Dor nebula. But in the SIS map, scattered X-rays can be considerable near N157B at energies 
^ 1.5 keV. 
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Fig. 3 presents close-ups of N157B in X-ray, compared with images in Ha and in radio 
continuum. The northeastern half of the PSPC image (Fig. 3a) shows a clear excess of diffuse X-ray 
radiation, but this diffuse enhancement is harder to see in the raw distribution of RHRI counts 
(Fig. 3b), which include a large number of non-cosmic X-ray events. The low surface brightness 
diffuse emission is, however, evident in the smoothed RHRI maps (Figs. 3c-d). The distinct part 
of the enhancement, outlined by the lowest solid contour level (7.4 x lO^^counts arcmin""^) in 
Fig. 3c, is somewhat square-shaped, and has a dimension of ~ 85" x 85". 

On top of the low surface brightness diffuse X-ray emission is a comet-shaped X-ray feature. 
The extended morphology of this feature, similar in the two RHRI observations, cannot be due 
to any instrumental effects. Point-like objects in the field do not show such elongation (Fig. 2; 
Wang 1995). This feature has a position angle of ~ 38° west to the north. At the southeastern 
end of the feature is a compact source {R.A. = 5''37™47^6;L>ed. = -69°10'20"), which appears 
extended in rh600228. A maximum likelihood fit to the count distribution of the source suggests 
count rate of (25 ± 5) x 10~^ counts s~^ and a PSF-subtracted 90% upper limit to the source's full 
Gaussian size as ~ 7". But an extent of a few arcseconds could be caused by the pointing "jitter" 
of ROSAT . Thus, the source is only marginally resolved, at the best. 

A different perspective of the data is presented in Fig. 4, where individual cuts of the RHRI 
intensity distribution are plotted. The compact source of the comet-shaped X-ray feature appears 
broader than the RHRI PSF (Fig. 4a), although part of this broadening may be due to the 
pointing uncertainty in the observation rh600228. The broad shoulder on the right (northwestern) 
side of the peak gives an estimate of the length of the X-ray feature as ^ 30" ± 6". The total 
width of the feature is ~ 20" ± 4" (Fig. 4b). Table 2 summarizes the properties of the individual 
X-ray components. 

3.3. Spectral Properties 

Figs. 5 and 6 show a remarkable similarity between the SIS spectra of SNR 0540-69.3 and 
N157B. Both spectra are flat and contain no apparent line features above 1.5 keV, setting these 
two remnants apart from young shell-like LMC SNRs studied by Hughes et al. (1995). In the 
SIS spectra of the shell-like SNRs, strong lines dominate, and weak continua fall steeply above 
~ 1 keV. 

The spectra of SNR 0540-69.3 and N157B show low energy turnovers at ~ 1 keV, indicating 
that the absorption is mainly due to oxygen and neon. The interstellar abundances of these 
elements are measured typically within the range between 30-40% solar in the LMC (Russell 
Sz Dopita 1992; de Boer et al. 1985), and around 40-60% solar in the solar neighborhood 
(e.g., Meyer 1995). Because the H I column density of the Milky way in the direction is only 
3.2 ± 1.3 X 10^° cm^^ (Bessell 1991), the absorption in the LMC dominates. We assumed an 
average metal abundance of 40% solar for X-ray-absorbing gas (Morrison &; McCammon 1983) 
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toward both N157B and SNR 0540-69.3. A slightly different choice of the metal abundance would 

only affect the equivalent hydrogen column density Nh in a spectral fit, but would hardly affect 
the fitting of other spectral parameters. Table 3 summarizes spectral parameters from various fits 
to the spectral data of the two remnants. 

The spectrum of SNR 0540-69.3 is well fitted with a simple power law. The power law slope 
is in good agreement with that obtained by Finley et al. (1993), whose analysis is based on a 
PSPC spectrum and assumes a solar metal abundance. Switching the abundance in our fit to 
the solar value and fixing the slope to their best-fit value 1.0, we obtained almost the same Nh 
and 0.1-2.4 keV band flux as theirs. This comparison suggests no gross inconsistency between the 
ROSAT and ASCA data. 

A power law, however, does not fit the spectrum of N157B as satisfactorily as that of SNR 
0540-69.3. The power law fit of the N157B spectrum can be rejected at the ^ 99% confidence. 
The residuals of the fit (Fig. 6) exhibit significant systematic fluctuations on energy scales of 
~ 0.1 — 0.3 keV, indicating possible emission line contributions from an optically-thin thermal 
plasma. The feature at ~ 0.9 keV, however, appears a bit too narrow, and its astronomical nature 
may be questioned, although a stronger line of similar width at ~ 1 keV has been reported in the 
spectrum of the nearby X-ray binary system 4U1626-67 (Angelini et al. 1994). Other possible line 
features are at ~ 1.4 keV and 1.8 keV. 

To probe the significance of the line features and the thermal plasma component in general, 
we made a joint spectral analysis of the SIS and PSPC spectra of N157B. First, we fitted the 
two spectra with a power law, allowing the model normalizations to vary independently. A 
significant excess around 0.9 keV is apparent in the PSPC spectrum (Fig. 7), even greater than 
that in the SIS spectrum. An inclusion of a simple Gaussian line improves the fit significantly, 
reducing the x^/d-o.f. from 197/150 to 171/146; the excess is gone in the residuals for both 
the SIS and PSPC spectra (Fig. 8; Table 3). The line in the PSPC spectrum has a flux 3.5 
(2.4-5.5) xlO""^ photons s"! cm'^ at 0.9(0.89 - 0.92) keV. The flux in the SIS spectrum is, 
however, only 45(21-68)% of that in the PSPC spectrum, whereas the power law normalizations 
in the SIS and PSPC spectra are very close (within ~ 5%). This line flux discrepancy can be 
explained if the thermal component originates in the low surface brightness diffuse emission region 
(Table 2). Part of the diffuse emission is not accounted for in our extraction of the SIS spectra, 
which assumes the remnant to be point-like. 

We further tried the standard collisionally equilibrium plasma model (i.e., Raymond &; Smith 
1994, R&S hereafter). For simplicity, we discuss here only the results from the fits to the SIS 
spectrum. A single-temperature model fit is not satisfactory; the best-fit model is much steeper 
than the SIS spectrum at energies ^ 4 kcV. Such a thermal model is, however, perfectly acceptable 
for the Einstein SSS data, partly because of the instrument's narrower energy coverage ^ 4 keV. 
The absorption inferred from the SSS spectrum is also significantly smaller than the value obtained 
from the ASCA data. We suspect that this is caused by a considerable soft X-ray contribution 
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from the main body of 30 Dor within the SSS aperture of ~ 6' diameter. Therefore, the results 
from the SSS spectrum should be used with caution. A two-component model consisting of a 
power law and an R&S plasma improves our power law fit to the SIS spectrum only slightly 
(x^/d.o./. = 164/125); the thermal component, even with the abundance allowed to vary, does 
not adequately account for the apparent line features. Without accurate information on the 
remnant's environment, structure, and evolution, we are unable to use the limited statistics and 
line signatures of the data to explore the vast parameter space of more sophisticated models. 

The power-law contribution may be determined, however, if it is solely due to the comet- 
shaped feature, presumably a synchrotron nebula (§4.3). We used the RHRI count rate of 
0.054 counts s~^ (Table 2; with an estimated uncertainty of ~ 20%) to isolate the contribution of 
the feature below 2 keV. At higher energies, where the thermal contribution may be neglected, 
we used the SIS spectrum. The spectral parameters from a power fit to this SIS/RHRI combined 
spectral data are included in Table 3. Because of the exclusion of the diffuse component, the Nh 
value is greater than those from the previous fits. The luminosity of the power law component is 
Lx = 1.0 X 10^^ ergs s~^ in the 2-7 keV band, or = 5.6 x 10^^ ergs in the 0.2-4 keV band. 
The uncertainties of these luminosity estimates are ~ 30%. 

Fig. 9 presents the power law model, together with the SIS spectrum including data points 
below 2 keV. The residuals, evidently showing line features, do appear to be thermal in origin. A fit 
with a single temperature R&S model suggests a characteristic plasma temperature of ~ 0.7 keV. 
However, the fit {x^ /d.o.f. = 187/128) is not satisfactory; there is a significant excess above the 
model below 0.7 keV. A similar fit {x^ /d.o.f- = 21/26) to the residuals in the PSPC spectrum is 
good, and the best-fit temperature is 0.41(0.31 — 0.53) keV. Thus, the thermal model describes the 
residuals well only at the PSPC's resolution. The predicted RHRI count rate ~ 0.031 counts s~^ 
is consistent with the rate of the diffuse component within an estimated systematic error of up to 
~ 40%. 

The above isolation of the power law component is physically plausible, but not unique. 
Alternatively, the power law component might arise only in the compact X-ray source. A power 
law fit to this combination (Table 3), however, is not as satisfactory. In particular, the absorption 
appears much too large (by a factor of ~ 3) to be consistent with the measurement of the optical 
extinction toward N157B (§4.4). 

4. Discussion 

4.1. Qualitative Considerations 

The striking similarity between the spectra of N157B and SNR 0540-69.3, together with the 
centrally-peaked morphology of N157B in both radio and X-ray, leaves little doubt about N157B 
as a Crab-like SNR. In contrast to the highly polarized emission observed in the Crab and SNR 
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0540-69.3, however, no radio polarization has been detected in N157B (D94). 

The comet-shaped X-ray feature is evidently the most interesting emission component of 
N157B. Markwardt &; Ogelman (1995) have observed a morphologically similar feature near the 
Vela pulsar. They interpret this feature as the thermal emission from a cocoon around a jet from 
the pulsar. The pulsar's proper motion is found to be nearly perpendicular to the direction of 
the jet. Therefore, it is possible that the jet is sweeping up the ambient material, which may 
explain the apparent higher density in the cocoon than in the ambient medium. However, the 
X-ray luminosity of the N157B feature is more than three orders of magnitude greater than that 
(~ 10^^ ergs s~^) of the Vela jet. Furthermore, the ASCA spectrum shows that the emission, most 
of which arises in the comet-shaped feature, is predominantly nonthermal in origin (§3.3). 

4.2. Low Surface Brightness Diffuse X-ray Emission 

The diffuse emission component of N157B most likely arises in a diffuse plasma, which is 
responsible for the possible weak line features in the ASCA spectrum. We adopt a conversion of an 
RHRI count rate to an emission integral (EI — defined as / dV) as 2.7 x 10^*^ cm~^/ (counts s~^). 
This conversion applies to the model of Nh ~ 10^^ cm~^ and T ~ 0.7 keV, and is good within 
a factor of 2 as long as the temperature ^ 0.4 keV. Prom Table 2, we infer a count rate of the 
diffuse component as 0.019 counts s~^, corresponding to an EI of 5.0 x 10^^ cm~^ and an X-ray 
luminosity as 1 x 10^^ ergs s~^ in 0.1-2.4 keV band, or 0.7 x 10^^ ergs s~^ in the 0.5-2 keV 
band. Further assuming that the diffuse component arises in an approximately spherical space 
of diameter D = 20 pc, and using ^ as the effective spatial filling factor, we estimate the gas 
density n/j as ~ 0.6C^°-^ cm~^. The corresponding total gas mass is then ~ 56C°'^ Mq. If C ^ 0.1, 
the remnant expansion should then be in the Sedov phase (e.g.. Cox 1972). The total thermal 
energy in the X-ray-emitting plasma is ~ 3 x IQ^^(^^-^ ergs, a considerable fraction of the expected 
mechanical energy of an SN. Because this fraction is expected to be ^ 72% in the Sedov phase, 
the inferred SN energy is in a nominal range of a few times 10^*^ ergs. Finally, the Sedov age of 
the remnant is tsNR = 5 x 10^D20pc^o TkcV y^^' where -D20pc is the remnant diameter in units of 
20 pc and To.ykeV the diffuse plasma temperature in units of 0.7 keV. 

Both the above X-ray size and boundary agree with those inferred from kinematic studies 
of Ha-emitting gas in N157B (e.g., CKSL92). These studies have revealed shocked materials 
with velocity offsets up to ~ 10^ km s~^. However, since the velocity of Ha-emitting gas does 
not necessarily represent the expansion of the outer shock of the N157B remnant, the inferred 
characteristic age can, at most, be considered as an upper limit to the true age of the remnant. 
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4.3. Comet-Shaped X-ray Feature 

The predominant non-thermal nature of N157B leads us to explain the comet-shaped X-ray 
feature as a synchrotron nebula, powered by the putative pulsar embedded in N157B. This 
explanation, constrained by the observations, provides specific predictions about the energetics 
and dynamics as well as the spatial and spectral properties of the remnant. 

Prom the luminosity of the X-ray feature, = 5.6 x 10^^ ergs in the 0.2-4 keV band (§3.3), 
we may estimate the spin-down power E of the pulsar. A linear relation logL-j « 1.39 log£'— 16.6, 
obtained by Seward &; Wang (1988), fits reasonably well to the three confirmed Crab-like 
SNRs. Applying this relation to N157B, we obtain E 2 x IQt'^ ergs s~^. Following Seward 
&L Wang (1988) again, we further estimate the period and magnetic field strength of the pulsar 
as ~ (0.08 s){tzE2.s)~^^'^ and ^ (1 x 10^^ G)t^'^ E-^^^"^ , where ts is the remnant's age in units of 
10^ yrs. 

The steep X-ray spectrum provides clues about the energy losing mechanism of relativistic 
particles in the synchrotron nebula. The X-ray-emitting lifetime tx of an electron in the nebula 
can be expressed as 

~ (40 yrs)e-V2(i7_4)-3/2^ (1) 

where e is the characteristic synchrotron photon energy in units of keV, and iy_4 is the magnetic 
field in units of 10~^ G. This timescale is much shorter than the expected age of the remnant. 
Therefore, a steady pulsar spin-down luminosity is a good approximation to model the current 
X-ray emission from the nebula. We assume that the distribution of pulsar injected particles is 
approximately a power law with an index of 7, and that the velocity distribution is isotropic. The 
X-ray emission spectrum is then a power law with a = 7/2 if synchrotron radiation loss dominates, 
or a = (7 — l)/2 if the energy loss is predominately due to adiabatic expansion of pulsar wind 
materials (e.g.. Tucker 1977). In general, one expects a convex-shaped spectrum: The adiabatic 
loss dominates at low energies, whereas the synchrotron becomes important at high energies. Such 
a slope change from high to low energies is a general characteristic of Crab- like SNRs (e.g., Du 
Plessis et al. 1995). The steep spectrum observed in N157B thus suggests that the synchrotron 
loss dominates in N157B even in the relatively low energy range 0.5-7 keV. Prom the slope of ~ 1.5 
(§3.3), we obtain 7^3. 

The spatial properties of the comet-shaped X-ray feature have implications for the dynamics 
of relativistic electrons. Because of the short lifetime of the electrons, the size of a synchrotron 
nebula should typically be much smaller than that of the corresponding radio one if electrons 
diffuse through the nebula medium at the Alfven velocity. The observed size and comet-like 
shape of the X-ray feature suggest that X-ray-emitting particles are undergoing an outflow, 
which can travel much faster than diffusing electrons. The compact source, where energy is 
likely being deposited by the pulsar, is separated from the radio peak (§4.4), where most of the 
pulsar spin-down energy has supposedly been dumped. If the birth place of the pulsar is close to 
the position of the radio peak and has an age of ~ 5 x 10^ yrs (§4.2), the separation of ~ 12" 
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between the pulsar and the radio peak then imphes that the pulsar has a transverse velocity of 

~ 6 X 10^ km s^^, compared to an average pulsar velocity of ~ 4 x 10^ km s^^ (Lyne &: Lorimer 
1994). We use Vp^^ (in units of 10^ km s~^) to characterize the total velocity of the pulsar motion. 

Such a relatively fast moving young pulsar can produce a strongly elongated pulsar wind 
bubble (Fig. 10). This bubble consists of several distinct regions. First, one may expect a bow 
shock running ahead of the pulsar, while the corresponding reverse shock terminates the free 
pulsar wind. The size of the region enclosed by this bow shock can be small, but the pressure 
can be large, depending on the ram pressure. Enhanced synchrotron radiation can then make 
the region especially bright in X-ray. Second, accelerated by the pressure gradient between the 
bow shock and the ambient medium, the shocked wind materials, which are relativistically hot, 
naturally forms a supersonic jet inside a tunnel in the opposite direction of the pulsar's motion 
(Wang, Li, &: Begelman 1993). Third, the jet likely shoots through the bubble, and then induces 
another bow shock (or a series of oblique shocks) on the other side. Fourth, the terminated wind 
materials are forced back to the bubble, where most of the pulsar spin-down luminosity is dumped. 
Let us now examine this picture in a more quantitative way. 

The size of the bow shock region around the pulsar should be consistent with the observed 
upper limit. The characteristic reverse shock radius is defined by the ram-pressure balance 
between the pulsar wind and the pulsar motion, 

r, = (4 X 10-2 pc)[E3s/{navls)]'^^ (2) 

where E^s is the spin-down luminosity in units of 10^^ ergs s~^ and Ua is the ambient medium 
density, which may be of the same order of magnitude as nfi estimated in the diffuse X-ray emission 
(§4.2). The size and shape of the outer bow shock depend critically on the dynamical effect of 
magnetic fields, and no theoretical model has yet been developed to describe this phenomenon. 
We thus characterize the total volume of the region enclosed by the bow shock as = |7r(^rs)^, 
where the parameter ^ is to be constrained. Assuming that the energy is partitioned equally 
between particles and magnetic fields, that the field is randomly distributed, and that the pressure 
is balanced in the region, we have 

H « (0.8 mG)inavls)'/^ (3) 
and ^ 

We = ^V^, (4) 

where We is the total energy in particles. Since only about half of the observed X-rays from the 
feature originates in the compact X-ray source (Table 2), the time spent by the shocked pulsar 
wind materials in the bow shock region is probably comparable to the lifetime of an electron that 
emits 1 keV photons. For a steady jet, we expect 



We - t^E = (5 X 10^4 ergs)£;38(na^;2 3)-3/4 



(5) 
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Prom (l)-(5), we obtain 

The bow shock has a size ~ 2^rs ~ 1" E^^^ {uaVp ^)''^^^'^ , consistent with the upper Umit to the full 
size of the compact X-ray source 7"; Table 2). 

After being forced out from the bow shock, the wind materials can continue to radiate in 
X-ray (Fig. 10). The radiation from the relativistic jet itself (Wang, Li, &; Begelman 1993), 

however, may not be significant. Most of the X-rays from the comet-shaped feature outside the 
compact X-ray source are probably released after the jet is terminated and are from shocked wind 
materials inside the bubble, explaining why the width of the X-ray feature (Table 2) is comparable 
to that of the compact radio emission region (§4.4). Just like particles in the bow shock region, 
the shock materials in the wind bubble are still relativistically hot, and the emission is primarily 
synchrotron. 



4.4. Radio Emission 

We speculate that the lack of polarization in the observed radio emission of N157B is due to 
the presence of a giant H II region produced by the coexisting OB association LH99 (Schild &: 
Testor 1992). CKSL92 find that about 20% of the radio emission arises in the H II region, to be 
consistent with both the Ha-to-radio continuum ratio and the average extinction {Ay = 1.0 mag), 
inferred from Balmer decrement measurements over a 4^9 region surrounding N157B (Caplan &; 
Deharveng 1985). The rest of the radio emission is presumably nonthermal; part of the emission 
should be associated with the synchrotron nebula of the N157B remnant. The radio emission 
above the level of 28 mJy beam~^ in Fig. 3d has an extent similar to the X-ray feature, but the 
centroid of the emission feature is about 12" it 3" away from the compact X-ray source. Based on 
the relations Ay/E{B — V) = 3.1 and Nh/E(B — V) = 2.4 x 10^^ cm~^, which are appropriate 
for the interstellar gas in the 30 Dor region (Fitzpatrick 1986), we infer Nh = 7.7 x 10^^ cm~^. In 
comparison, our X-ray- measured Nh is ~ 1 x 10^^ cm~^ (§3-3). Thus the N157B remnant is likely 
located in or behind the H II region. 

The ionized gas and the expected magnetic field of the H II region can cause Faraday rotation 
in the radiation from the synchrotron nebula. Assuming that the 20% of the observed radio 
emission arises in a spherical region of ~ 20 pc diameter, similar to the projected dimension, we 
estimate a mean free electron density rig ~ 50 cm^'^. The Faraday rotation angle varies with 
ne-ff||, where ify is the line-of-sight intensity of the field. Following Burn (1966), we express the 
dispersion of the rotation within a radio survey beam as 

6(1) ~ [KneH^^deX^] ds/de, (7) 

where the coefficient K = 2.6 x 10~^^, de is the characteristic size of turbulent eddies, dg = 1.8 pc 
and A = 13 cm are the beam size and the wavelength of the 2.3 GHz survey. In the bracket of the 
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above equation is the Faraday rotation angle through a single eddy. Assuming a typical interstellar 
field of ~ 5/xG, we obtain 5(f) ~ 6. This dispersion produces a radiation depolarization 

p ^ Pi exp[-2{54ifdHn/de], (8) 

where de should be comparable to, or smaller than, the dimension dun of the region. Therefore, 
only if a fraction of the observed H II region is in front of the nebula can the Faraday dispersion 
explain the lack of the polarized emission in the 2.3 GHz survey. Clearly, a radio survey at a 
higher frequency will have a better chance to detect the polarized emission from the synchrotron 
nebula. 

5. Summary 

Based on the ROSAT and ASCA observations, wc have examined the X-ray timing, spatial, 
and spectral properties of N157B. An algorithm has been devised for extracting ASCA spectra 
for a point-like X-ray source. For comparison, we have also presented an ASCA spectrum of SNR 
0540-69.3, the only LMC SNR that contains a detected pulsar. We have arrived the following 
main results and conclusions. 

1. The fiat and featureless spectrum of N157B in the 2 — 7 keV range resembles that of SNR 
0540-69.3, and is distinctly different from those of young shell-like SNRs in the LMC. The 
power law energy slope (~ 1.5) of the spectrum is, however, steeper than those 1.0) of 
SNR 0540-69.3, MSH 15-52, and the Crab nebula. 

2. The ASCA spectrum of N157B shows possible line features at energies ^ 2 keV. A joint 
analysis of the ASCA and ROSAT data suggests that these features are likely associated 
with a thermal emission component of the remnant. 

3. No pulsed signal is identified; the pulsed fraction is <^ 17% of the remnant's total luminosity 
in the 0.5-2 keV band, or ^ 9% in the 2-10 keV band. So the putative pulsar in N157B is 
still elusive. 

4. More than two thirds of the observed X-ray radiation from N157B comes from a comet-shaped 
X-ray feature, which has a size of '-^ 7 pc x 5 pc and a luminosity of ~ 1 x 10^^ ergs s~^ in 
the 2-7 keV range. 

5. We interpret the X-ray feature as the synchrotron radiation from pulsar wind materials. 
The pulsar is likely moving at a speed of ^ 10'^ km s~^. A ram-pressure confined bow shock 
around the pulsar may account for an embedded compact X-ray source, which is probably 
unresolved. The ram pressure further focuses fresh wind materials into a relativistic jet. The 
jet is eventually terminated on the other side of the pulsar wind bubble, where most of the 
spin-down energy has been dumped. 
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6. Around the X-ray feature, we further detect a low surface brightness diffuse X-ray emission 
region. This region most hkely accounts for the apparent thermal component in the 
ROSAT and ASCA spectra below ~ 2 keV. The remnant has a diameter of ~ 20 pc and a 
0.5-2 keV luminosity of ~ 7 x 10^^ ergs s~^, a total thermal energy of ~ 3 x lO^'' ergs, and 
an age of ~ 5 x 10^ yrs. 

7. The lack of polarized radio emission from N157B may be explained by the depolarization of 
a foreground H II region. 

We thank J. Dickel for the radio map, C. Smith for the Ha image used in the paper, and 
G. Testor for sending us his optical CCD images of the N157B region. We thank the referees 
for various critical comments that led to improvements in the paper. This work is funded by a 
Lindheimer Fellowship and by NASA grant NAG5-2717. 

A. ASCA Background Spectral Estimation for an Isolated Point Source 

A.l. Introduction 

Estimating a local background spectrum and subtracting it from an on-source X-ray spectrum 
normally represents a crucial step in analyzing the spectral properties of a source. Because of 
the broad ASCA point spread function (PSF) (sec §2.1), it is not usually practical to obtain an 
independent measurement of the background free of source contamination without going far from 
the source, where data might not exist, or the background might differ significantly. For example, 
an annular background region of radius 3' - 5' would require that the source count rate be less 
than 0.03 counts per second for it not to contribute more than 10% to the total counts in the 
background region. This count rate is barely detectable in a typical 20 ks observation. For an 
observation of a moderately bright point source, the counts in the background are an admixture 
of source and background counts. On the other extreme, a source of a few count per sec would 
dominate the background region, contributing > 90%, and the background can be ignored to 
first-order. We are interested here in the intermediate case. 

The measured background includes X-rays of cosmic origin. X-rays and particles of Solar 
origin, cosmic ray induced X-rays, and internal instrumental counts. Perhaps most insidious for 
ASCA is the stray light X-rays (mirror rays which undergoes single, instead of double reflections, 
etc.), which are not imaged in the usual way, and which depends strongly on the celestial 
distribution of X-ray emission in and out 1°) of the telescope's field of view. The background 
depends on the satellite geographical location and height above the Earth, solar activity, and 
the incident cosmic X-rays as well as telescope pointing angle relative to the Earth, Sun, and 
Earth's magnetic field vector. Taken together, these produce a background for a given observation 
whose spectrum and intensity depends on date, orbit, sky, and detector coordinates. The broad 
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PSF wings also mean that other sources in the field may substantially color the target spectra, 
producing an additional effective background. 

The current recommended practice for estimating the background is to extract counts from 
the ASCA Deep Field pointings, screened in a manner identical to the data, and selected from 
the source detector region. For the simple case of an isolated point source, using the Deep 
Fields should give a reasonable zero-th order estimation of the background. However, for a given 
observation and target, this practice is limited by the uniqueness of background as discussed above. 
As an added complication, although these fields might represent the background statistically in an 
average field, they are known to contain explicit point sources, which may compromise extracted 
spectra near their detector locations. Furthermore, the response of the detector to the background 
changes over the course of the mission due to known secular degradation of the instruments 
with time. Most importantly, the background intensity and spectrum can vary from one field to 
another, especially in regions of nearby galaxies such as the LMC. 

Here, we present a simple method for a first-order estimate of both the source spectrum and 
the background spectrum from the data itself, for the case of an isolated point ASCA source 
embedded in a statistically uniform local background. The advantage of self-calibrating is evident 
given the strong temporal, instrumental, orbital, and celestial dependence of the background. Our 
method is to consider the spectral ratio of observed and expected counts in concentric annuli 
centered on the source, each of which is an admixture of the true source and background spectrum. 
We look for deviations of the radial average profile from that expected for a point source and 
simultaneously solve for the two spatially distinct spectral components. Thus we can estimate the 
background over the source region itself. This is in effect a form of spectro-spatial deconvolution 
of the background and source spectra. 

A.2. Method 

Consider the observed counts in two concentric regions centered on the source and in the 
measured spectral bin ("PHA" bin) E' , Oi{E') and 02{E'), both of which are a composition of 
the source and background spectra, S{E') and B{E') respectively, 

Oi{E') = pi{E)S{E') + aiB{E') (Al) 

and, 

02(^0 = P2{E)S{E') + a2B{E')v{E) (A2) 

where ai, are the geometric areas of the respective regions; v{E) is the mean vignetting of 
the background region relative to the source; p\{E) and P2{E) are the encircled energy function^ 



^The unfortunate use of the term "energy" here is not to be confused with the spectral energy, as it in fact refers 
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v(EEFs) within those areas. Here, the energy, E, refers to the energy dependence of p{E) 
and v{E), which also depend on specific locations of ai and 02 in detector coordinates. These 
dependencies have been modeled with laboratory and on-orbit calibrations (see §A1.3, below) and 
therefore, the ratios 7 = ai/v{E)a2, and 5 = pi/p2, are determined. Notice that these ratio is not 
strongly energy dependent below ~ 6 keV, and for most cases, to the limit of our approximation, 
a single number for each is reasonable. 

We can then solve for the source and background spectra, and find. 



where piS{E') is the total counts enclosed in the source area ai in spectral bin E' , and aiB{E') 
is the background spectra over the same area at that energy. The coefficients of S{E') and B(E') 
insure correct normalization and units. The resulting errors for the source and background are 
simply found from the quadratic sum of the errors associated with individual terms. 

In practice, we carry out this calculation on each PHA (PI) bin in our spectral file. The energy 
dependencies (E) and {E') are related by the spectral redistribution file, RMF(ii^;£''), and do not 
enjoy a one-to-one correspondence. However, our method is insensitive to this redistribution, as 
it is a third-order effect because it arises in the errors in the ratios 7 and 5, which come from the 
fraction of redistributed counts at the lower energies. The typical negative sloped spectra and the 
weak energy dependence of these ratios insure that this effect is negligible for our purpose. 

For simplicity, we assume the ratios 7 and 6 are energy independent. We can then easily 
compute the source and background spectral components using the FTOOL MATHPHA to 
perform the calculation on each element of the spectral ( "PHA" ) file. We are then left with a 
spectral file which can be imported directly into a spectral analysis package such as XSPEC. 
These spectra will be correctly normalized and the instrument and mirror response files ("RMF" 
-|- "ARF" )0 created for Oi [E') are still appropriate. It is most practical to use the observed source 
spectra Oi{E') along with the derived background spectra aiBi[E') in analyzing the source 
spectra. 

The energy dependence of p{E) and v{E) can be explicitly accounted for by using the Awaki's 
"XRTEA" effective area subroutine which returns p{E) and v{E) as a function of source position 
in a detector and off-source angle. This code is available in the FTOOLS software package. An 
attractive alternative for measuring the EEF is to use an observation of a bright point-like source. 
If this fiducial source is acquired at the same detector position (i.e. 1 CCD mode default) as 



to the power contained in the mirror PSF at various radii; it is defined as the relative counts ("energy") enclosed in 
the PSF between the radii a and b, EEF{E; a, b) = £ PSF{E) / PSF{E). 

description of the "PHA", "ARF", & "RMF" files, the XSPEC spectral analysis program, and the FTOOLS 
software package is available in the ASCA data analysis guide (Day et al. 1995) and references therein. 



Pi{E)S{E') = 5{Oi{E') - j02{E'))/{5 - 7) 



(A3) 



aiBiE') = -i{Oi{E') - 502iE'))/{j - 6) 



(A4) 
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the source under study, and the data screened and exposure-corrected in the same way, then 
the backgroTind for the fiducial source relative to the study source may be reasonably ignored. 
Modeling the EEF with a real observation allows for the jitter in the attitude, which is similar 
from observation to observation (order ^ 10"). For the SIS, the fiducial source should not be too 
bright as to be noticeably affected by pile-up. The bright CVs make excellent targets, and have 
been successfully used for this purpose. 

In conclusion, the above method allows an estimate of the background spectrum under a 
circular aperture, centered on a point-like source. Given enough counts, this method can naturally 
extend to a least square solution for a general spectro-spatial deconvolution for a point source. We 
may extract counts from a set of N concentric annuli and solve the set of linear equations of the 
form: 

On{E') = Pn{E)S{E') + anB{E')vn{E) (A5) 



Table 1: Summary of ASCA Exposures on N157B 



Detector 


Exposure 


Source 


Off-axis 


Vignetting 


(CCD Chip) 


Time 


Counts'^ 


Distance 


Fraction 




(s) 




(arcmin) 






ASCA Seq. 


No. 20000000 




SISO (C2) 


16743 


2961 


3.06 


0.90 


SISl (CO) 


19245 


2717 


7.17 


0.64 


GIS2 


30820 


6784 


3.43 


0.87 


GIS3 


30816 


5856 


7.14 


0.64 




ASCA Seq. 


No. 50002000 




SISO (C0+C3) 


26123 


2758 


8.97 


0.54 


SISl (C2+C1) 


16876 


2623 


4.01 


0.84 


GIS2 


30543 


5727 


7.24 


0.64 


GIS3 


30538 


6171 


7.12 


0.64 



" SIS and GIS counts are obtained within I'.l and 4!0 radius circles centered on N157B. 



Table 2: Spatial Components of N157B 



Component 


X-ray Full Size 


Count Rate*^ 


Comment 


Whole Remnant 


85" X 85" 


73 


^ 7.4 X 10""^ counts s~^ arcmin"^ 


Comet-Shaped Feature 


30" X 20" 


54 


^ 2.2 X 10~^counts s~^ arcmin"^ 


Compact source 


^7" 


25 


Probably unresolved 



"in units of 10 ^ counts s 
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Table 3: Model Fit to the ASCA and ROSAT Spectral data'' 



Data 



Model 



'^/d.o.f. Energy slope Column Density 



Norm'' 



SNR 0540-69.3 



SIS 



Power law 



172/160 1.04(0.98-1.10) 7.9(7.4-8.4) 



11(10-12) 



N157B 



SIS-I-PSPC Power law+line 

SIS-(-RHRI(feature) Power law 
SIS-l-RHRI(source) Power law 



171/146 1.4(1.3-1.5) 6.9(6.3-7.5) 
82/78 1.5(1.4-1.7) 11(7.3-16) 
92/78 1.9(1.7-2.1) 28(22-35) 



3.2(3.0-3.4)^= 

3.5(2.8-4.5) 

6.1(4.8-7.4) 



"The parameter limits are all at the 90% confidence level. 
'At 1 keV and in units of 10"^ photons s"^ keV"^ cm"^. 
'^From the normalization of the PSPC spectrum. 
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Fig. 1. — ASCA SIS maps of the 30 Dor region. The maps are corrected for exposure and vignetting. 
The contours are at (3.0, 4.0, 5.5, 7.8, 11, 16, 24, 35, 52, 77, and 116) xlO"^ counts s"^ arcmin-^ 
in the 0.4-1.5 keV band, and at 1.7, 2.3, 3.2, 4.6, 6.6, 9.6, 14, 21, 31, 47, 69, 104, and 156) xlQ-^^ 
counts arcmin"^ in the 1.5 — 7 keV band. Faint radial strips near the southeastern boundaries 
of the map, especially in the 1.5 — 7 keV band, are due to stray light (reflected mirror rays) from 
LMC X-1 to the south. 

Fig. 2. — N157B and its environment: RHRI X-ray intensity contours overlaid on an Ha image. 
The X-ray intensity is calculated by coadding the two background-subtracted RHRI observations 
and is adaptively smoothed with a Gaussian of adjustable size to achieve a constant local count- 
to-noise ratio of ~ 8 over the image. The X-ray contours are at (2.1, 3.8, 5.9, 8.4, 12, 16, 22, 30, 
60, 120, 250, 450, and 800) xlO"^ counts s^^ arcmin-^. 

Fig. 3. — Close-up views of the N157B region: the PSPC count distribution in the 0.5-2 keV band 
(a), the RHRI count distribution (b), the RHRI intensity contours compared with an Ha gray- 
scale map (c), and 2.3 GHz continuum intensity contours overlaid on the gray-scaled RHRI map 
(d) . In (c) , the X-ray intensity background-subtracted and exposure-corrected, is smoothed with a 
count-to-noise ratio of ~ 4, and the contours are at (3.6, 7.4, 13, 22, 35, 54, 120, 250, 450, and 800) 
xlO^^counts arcmin^^. In (d), the gray-scaled X-ray intensity is plotted logarithmically in the 
range between 2.4 x 10~^ (about la above the local background) and 0.5 counts arcmin"^, and 
the 2.3 GHz radio continuum contours are at 2, 4, 8, 16, 28, 44, 64, and 88 mJy beam~^. The radio 
data (see also D94) have a FWHM beam size of ~ 7'.'5, comparable to that of the RHRI image. 

Fig. 4. — RHRI X-ray intensity distribution in strips across the field of N157B. Each strip is 15" 
wide. The histograms are calculated so that each bin has at least 25 counts. Fig. 4a presents 
a single SE-NW cut through the feature (38° west to the north), together with the RHRI PSF 
(the Gaussian- like profile) and the radio intensity distribution (solid curve) in the same cut; the 
southeast is to the left. The label of the ordinate is for the X-ray data, whereas the radio intensity 
is in units of Jy beam~^). Fig. 4b includes the NE-SW oriented strips, perpendicular to the N157B 
feature. The northeast is to the left, and the most northwestern strip is at the top. The fourth 
strip is centered on the position of the compact X-ray source. The bottom histogram is in units as 
posted; for ease of separation, the intensity in the subsequent histogram is multiplied by 10, and 
so on. 

Fig. 5. — ASCA SIS spectrum of SNR 0540-69.3. The best-fit power law model is presented as the 
histogram. The residuals (spectrum minus model) are shown in the lower panel. 

Fig. 6. — ASCA SIS spectrum of N157B. The rest is the same as Fig. 5. 

Fig. 7. — Joint power law fit to the ASCA and PSPC spectra of N157B. The rest is the same as 
Fig. 5. 



- 21 - 



Fig. 8. — Joint fit to the ASCA and PSPC spectra of N157B with a two-component model consisting 
of a power law and a Gaussian line. The rest is the same as Fig. 5. 



Fig. 9. — ASCA spectrum of N157B and the best joint-fit of a power law to both the spectrum 
above 2 keV and the RHRI flux of the comet-shaped X-ray feature (§3.3). The RHRI data point 
has its energy coverage between the 0.3-2.4 keV range. The residuals of the fit are shown in the 
lower panel. 



Fig. 10. — Illustration of the current pulsar position and the flow of pulsar wind materials 
interacting with the N157B remnant. The arrows represent the flow of the materials. The remnant's 
blast wave is represented by the outer shell. The configuration of the components is somewhat 
similar to the fan-shaped radio feature connecting the pulsar PSR 1757-24 and the SNR G5.4-1.2 
(Frail & Kulkarni 1991). 
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